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Abstract 

The inhomogeneous distribution of matter in the non-linear regime of 
galaxies, clusters of galaxies and voids is described by an exact, spherically 
symmetric inhomogeneous solution of Einstein's gravitational field equations, 
corresponding to an under-dense void. The solution becomes the homoge- 
neous and isotropic Einstein-de Sitter solution for a red shift z > 10 — 20, 
which describes the matter dominated CMB data with small inhomogeneities 
5p/p ~ 10~ 5 . A spatial volume averaging of physical quantities is introduced 
and the averaged time evolution expansion parameter 9 in the Raychaudhuri 
equation can give rise in the late-time universe to a volume averaged decel- 
eration parameter (q) that is negative for a positive matter density. This 
allows for a region of accelerated expansion which does not require a neg- 
ative pressure dark energy or a cosmological constant. A negative decel- 
eration parameter can be derived by this volume averaging procedure from 
the Lemaitre-Tolman-Bondi open void solution, which describes the late-time 
non-linear regime associated with galaxies and under-dense voids and solves 
the "coincidence" problem. 
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1 Introduction 

In a recent article, we investigated a cosmology in which a spherically symmetric in- 
homogeneous enhancement is embedded in an asymptotic FLRW universe jTj. The 
inhomogeneous enhancement is described by an exact inhomogeneous solution of 
Einstein's field equations. We found that the inhomogeneities can lead to a reinter- 
pretation of the luminosity distance di of a cosmological source in terms of its red 
shift z, owing to the observer dependence of these quantities. The time evolution and 
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the expansion rate of the inhomogeneous distribution of matter can lead to intrinsic 
effects such as cosmic variance at large angles and long-wavelength perturbations not 
described by a Friedmann-Lemaitre-Robertson- Walker (FLRW) homogeneous and 
isotropic universe. Therefore, the interpretation of the data using a FLRW model 
that the accelerating expansion of the universe is caused by dark energy may be mis- 
leading. This is important, for it is difficult to explain theoretically the postulated 
dark energy that causes the acceleration of the universe. The model also leads to an 
axis pointing towards the center of the spherically symmetric large scale inhomoge- 
neous enhancement with dipole, quadrupole and octopole moments aligned with the 
axis. It was shown that the luminosity distances and red shifts observed by different 
observers located at spatially different points of causally disconnected parts of the 
universe can have varying values. A spatial average of all these observations leads to 
an intrinsic cosmic variance in e.g. the deceleration parameter q. The distribution 
of CMB temperature fluctuations can be unevenly distributed in the northern and 
southern hemispheres. 

The acceleration of the expansion of the universe deduced from Type SNIa super- 
nova observations and the CMB WMAP data El Ej has been interpreted as due to 
the cosmological constant (vacuum energy), modifications of Einstein's gravitational 
field equations at large distances [Sj, and quintessence fields JB]. The quintessence 
explanations postulate a new form of matter with negative pressure called dark en- 
ergy. Recently, it has been suggested that the acceleration is caused by very long 
wavelength, super-horizon perturbations generated by a period of inflation in the 
early universe El El UD] ■ The back-reaction of perturbations on a FLRW back- 
ground universe has been the subject of investigation by several authors [H]. The 
predictions based on perturbation theory are limited by the condition $ 1, where 
$ is the gravitational potential and the assumption of a homogeneous and isotropic 
FLRW background. Moreover, the perturbation calculations do not possess a rig- 
orous gauge invariance unless all physical quantities are scalars under coordinate 
transformations. Otherwise, calculations rely on using specific gauges such as the 
synchronous gauge or the Poisson (Newtonian) gauge, which can possess residual 
gauge ambiguities. 

The suggestion that the acceleration of the universe can be caused by sub-horizon 
or super-horizon modes due to a phase of early universe inflation has been criti- 
cized (121 EH EU IIHI EE HZI ■ It is pointed out that second order perturbation effects 
of the form $V 2 $ are described by a renormalization of the local spatial curvature 
and cannot (for a positive energy density and pressure) produce a negative decelera- 
tion parameter. Hitara and Seljak ^3] use a derivation of the deceleration parameter 
q obtained from the Raychaudhuri ^Hl EH] equation that requires significant con- 
tributions from vorticity to produce a negative q when the strong energy condition 
p + 3p > is satisfied. Notari |Hj argues that the perturbation expansion fails at late 
times and leads to significant infrared effects and a possible negative deceleration 
parameter q. All of these arguments depend on a perturbation calculation against a 
homogeneous and isotropic FLRW background. Thus, the results are inferred from 
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the strong assumption that the universe is approximately homogeneous at all scales. 
We know that for distances of order < 150 Mpc there are galaxies, clusters of galax- 
ies and large scale voids ~ 50 Mpc, and that Sp/p > 1, so that linear perturbation 
theory fails. 

The inhomogeneous distribution of matter at late times could possibly allow for 
a negative deceleration parameter q, depending on the nature of the exact inhomo- 
geneous solution of Einstein's field equations. This will allow for the possibility of 
explaining the acceleration of the universe, without a negative pressure dark energy 
fluid or a cosmological constant. It can solve the "coincidence" problem, because 
the negative deceleration parameter q and the accelerating expansion will not occur 
in the absence of large-scale structure in the form of galaxies, voids and clusters of 
galaxies 1 

We shall describe the inhomogeneous late-time universe by an exact spherically 
symmetric solution of Einstein's field equations. The solution becomes the homo- 
geneous and isotropic Einstein-de Sitter matter dominated solution for z > 20, so 
that it can describe the WMAP data near the surface of last scattering with small 
temperature perturbations ST/T ~ 10~ 5 . However, for z < 20 the solution is inho- 
mogeneous and for luminosity distances corresponding to z ~ 0.2 — 1.4, includes the 
Type SNIa supernova measurements. We investigate the expansion geometry of the 
spherical distribution of inhomogeneous matter and under-dense voids [20] to see 
whether it allows for the possibility of a negative q and an accelerating expansion 
of the universe, without negative pressure dark energy or a cosmological constant. 

A generic relativistic solution of a three-dimensional inhomogeneous universe 
is presently unknown. The highly symmetric Lemaitre-Tolman-Bondi (LTB) [231 
I2U 122] solution that we shall consider in the following is inhomogeneous in only one 
spatial dimension. Even though the assumption of spherical symmetry is unrealistic, 
the model is useful in illustrating important local features of physical quantities. 

We will discover that although the LTB model appears to allow for a negative 
deceleration parameter for a positive density p > 0, an investigation of the gauge 
invariant Raychoudhuri equation that determines the rate of change of the volume 
expansion shows that for the synchronous and co moving LTB metric, the vorticity 
uj vanishes and q must be positive when the strong energy condition is satisfied. 
However, for our inhomogeneous late-time universe model, we are required to carry 
out a volume averaging of physical scalar quantities, such as the expansion parameter 
9. We find for such an averaging process that for a zero cosmological constant and 
for irrotational matter the average deceleration parameter (q) can be negative and 
describe an accelerating local region of the universe, corresponding to our observed 
Hubble radius, without invoking dark energy with a negative pressure. 

1 The possibility that a large-scale inhomogeneity can mimic a cosmological constant has been 
considered previously by Celerier |21j . The possibility that an under-dense void can cause the 
acceleration of the universe has been proposed by Tomita |22| . 
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2 Inhomogeneous Friedmann Equations 

For the sake of notational clarity, we write the FLRW line element 

dr 2 



ds * = d e - a 2 (t) ( + r 2 dtf) . 

\1 — kr 1 J 



where k = +1,0,-1 for a closed, flat and open universe, respectively, and dfl 2 = 
d8 2 + sm6 2 dcf) 2 . The general, spherically symmetric inhomogeneous line element is 

given by [2312312312312111211123111: 

ds 2 = dt 2 - X 2 (r, t)dr 2 - R 2 (r, t)dtt 2 . (2) 

The energy-momentum tensor T^ v takes the barytropic form 

T» v = {p + p)u»u v -p8» Vi (3) 

where m m = dx^/ds and, in general, the density p = p(r,t) and the pressure p = 
p(r,t) depend on both r and t. We have for comoving coordinates u° = l,u % = 
0, (i = 1, 2, 3) and g^ v u^u v = 1. 

The Einstein gravitational equations are 

+ = -SvrGTV, (4) 

where G^ u = R^ u — ^g^TZ, 1Z = g^ u R^ u and A is the cosmo logical constant. Solving 
the Gq\ = equation for the metric (j2J), we find that 

X(r,t) = ^M, (5) 
f( r ) 

where R' = dR/dr and f(r) is an arbitrary function of r. 
We obtain the two generalized Friedmann equations pQ: 

R 1R 2 11. , 2 , IR'R Iff 47rG. „ , 1. 
^ + 3^ + 3^ (1 - /} -3^^ + 3^ = -^ (P + 3p) + 3 A ' (?) 

where i? = dR/dt. 



3 Late-Time Matter Dominated Universe 

The late-time matter dominated universe will be pictured as a large-scale inho- 
mogeneous enhancement that is described by an exact inhomogeneous spherically 
symmetric solution of Einstein's field equations. The inhomogeneous enhancement 
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is embedded in a matter dominated universe that approaches asymptotically an 
Einstein-de Sitter universe as t — ► oo. An observer will be off-center from the 
origin of coordinates of the spherical inhomogeneous enhancement. The local inho- 
mogeneous solution corresponds to an open, hyperbolic solution and describes an 
under-dense void. Surveys such as the 2-degree Field Galaxy Redshift Survey and 
the Sloan Digital Sky Survey show large volumes of relatively empty space, or voids, 
in the distribution of galaxies [20J. 

For the matter dominated Lemaitre-Tolman-Bondi (LTB) j2HJ EH EHj model 
with zero pressure p = and zero cosmological constant A = 0, the Einstein field 
equations demand that R(r, t) satisfies 

2R(r,t)R\r,t) + 2R(r,t)(l - / 2 (r)) = F(r), (8) 

with F being an arbitrary function of r of class C 2 . 
The proper density of matter can be expressed as 

F' 

P = WnGR'R 2 ' ^ 

We can solve (jHJ) to obtain 

n-i s ^-i = 4-fiz/!_ 2 /A (io) 

p c H 2 S \ R 2 RR'J' 1 ; 

where 

H 2 S = H 2 ± + 2H L H r , (11) 
is an effective Hubble parameter and 

* = flx = i (12) 
We have for the critical density 

8nG Pc = H 2 S . (13) 

There are three possibilities for the curvature of spacetime: 1) f 2 > 1 open (fl — 1 < 
0), 2) f 2 = 1 flat (Q - 1 = 0), f 2 < 1 closed (O - 1 > 0). 

We describe the inhomogeneous density regime by a hyperbolic, open solution 
}' 2 {r) > 1 with p = A = 0, corresponding to an under-dense void. We choose 



f{r) = VI + r 2 (14) 

and the metric reduces to 

rfs 2 = dt 2 _ ^!M dr 2 _ ^r ^. (15) 
1 + r l 
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A parametric solution is given by 



R(r,t) = ^F(r)(f 2 (r) - l)' 1 (cosh u(r,t) - 1), (16) 
t + P(r) = ^(f(r) - l)- 3/2 (smhu(r,t) - u(r,t)), (17) 



8%Gp(r, t) 



(2F'(r)\( f 2 (r) - 1 



\F 2 (r)J\ R'(r,t) J Vsinh 4 \u(r, t). 

Here, (3(r) is an arbitrary function of r of class C 2 [25] , 

The function (5{r) can be specified in terms of a density on some spacelike hy- 
per surface t = to- The metric and density are singular for 



u 



or R'(r,t) = 0, (19) 



where u = refers to the hypersurface t + /3(r) = and the singular surface for 
R'(r,t) = is more complicated. Our pressureless model requires that the singular 
surface t(r) = E(r) describes the surface on which the universe becomes matter 
dominated (in the FLRW model this occurs at z ~ 10 4 ). The hyperbolic FLRW 
model is obtained by choosing the conditions 

f(r) = l + r 2 , F(r) = Asr 3 , p(r) = 0, (20) 

where s is a positive constant. This leads to the metric solution 

R = sr(coshu - 1), (21) 

t = s(sinhu — u), (22) 

8irGp= -^f , * V (23) 
4s 2 Vsmh 6 iM/ 

We shall change the radial coordinate to 

r* = r - r 5L 5, (24) 

where r = vsls is the location of the surface of last scattering corresponding to 
r* = at a red shift z ~ 1100. We postulate that (3(r*) approaches a maximum 
constant value at the approximate end of galaxy and cluster formation, Zf ~ 1 — 5, 
where z = Zf denotes the final value of the red shift. We assume that (3(r*) — > 
and fir*) — > 1 as r* — > 0. A model for /? is of the form: 



^)=^j. ( 25 ) 

where a, b and a are positive constants, (3(0) = and (3(r*) — ► a for r* ^> 6. The 
surface r* = for the metric is a singular surface. 
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The metric for f(r*) — > 1 becomes |2| 

d s 2 = dt 2 - (t + (3(r*))^ 3 (Y 2 (r*, t)dr* 2 + r*W), (26) 



where 



and 



y ™ = 1 + wWty (27) 

P ^ = 67rG(t + /3(r*)) 2 F(r*,t)' ^ 

We obtain for /3(r*) — > the homogeneous and isotropic solution corresponding to 
a spatially flat Einstein-de Sitter universe with the metric 

ds 2 = dt 2 - a 2 (t)(dr* 2 + r*W), (29) 

and a(t) oc t 2 ^ 3 . This solution is compatible as a background spacetime with the 
approximately homogeneous and isotropic WMAP data for 20 < z < 10 3 and Sp/ p ~ 

10- 5 m. 

The Einstein-de Sitter solution must be supplemented by non-zero pressure con- 
tributions for z ~ 10 3 in the linear regime at the surface of last scattering in or- 
der to take into account the radiation density pr. Moreover, a more general solu- 
tion is required that has a non-vanishing cosmological constant A, so that we have 
Q = Qm + + = 1 in agreement with the spatially flat result Q = 1.02 ± 0.02, 
obtained from the determination of the first acoustic peak in the WMAP power 
spectrum However, in the following section, we shall replace Q\ by a contribu- 
tion corresponding to an acceleration parameter q obtained from our inhomogeneous 
void solution with A = 0. 



4 The Inhomogeneous Late-Time Density and De- 
celeration Parameters 



We obtain from (JHJ) for f(r) = yl + r 2 and p = (we replace for convenience r* by 
r): 

Hi + 2H r H ± - — - — = SrrGp + A. (30) 
Dividing this equation by H\ gives 



ixGp A r 1 2r 2H, 



n ~ H 2 + R 2 + Rm 2 + R'RH 2 H\ ~ L (31) 



Let us expand R(r, t) in a Taylor series 
R(r, t) = R[r, t - (t - t)] = R(r, t ) 



l-(t -t)^^ + kto-t) 2 ^ to) 



R(r,t ) 2 V ' R(r,t. 



o 
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R{r,t c 



1 - (to - t)H ± - -{to - t) 2 q(r, to)H% x - 



(32) 



where to denotes the present epoch and H ± = R(r,to)/R(r,to). Moreover, we have 

For the matter dominated era we set p = in Eqs. © and (J7J) and substitute for R 
from (J7J) to obtain 

( + \ 1 i 47r PoMo) 
^^^SH 2 ^) 

1 r A l H 0r {r,to) 1 r 2 

3i?M \ 3^ ± (r,t ) 3H 0± {r,t ) 3 i? 2 # Y 1 j 

where #or(r,* ) = R'(r,t )/R'(r,to). 

If we set H ±(r, t ) = H 0r (r, t ) = H(t ) where H(t ) = a(t )/a(t ) and R(r, t ) = 
a(to), then we obtain the global FLRW expression for the deceleration parameter for 
a spatially flat universe with f 2 (r) = 1: 

<?0 = ~^0Af — ^oa, (35) 

where from (|3~T]) fioM + ^oa = 1 with Oom = ^Pom/^H 2 ^) and ^oa = A/3i/ 2 (io)- 
The variance of g is given by the exact non-perturbative expression: 

var(g) = (g 2_ (g) 2 ) l/2 = - (36) 

where 



( -» = "t^P (37) 

denotes the ensemble average and 7 denotes the determinant of the 3-dimensional 
metric gij (i = 1, 2, 3). 

We see from (I34J) that different observers located in different causally discon- 
nected parts of the sky will observe different values for the deceleration parameter 
q, depending upon their location and distance from the center of the spherically sym- 
metric inhomogeneous enhancement. This can lead to one form of cosmic variance, 
because the spatial average of all the observed values of local physical quantities, 
including the deceleration parameter, q, will have an intrinsic uncertainty. 

Let us set the cosmological constant to zero, A = 2 . From (pTTf and (jBH), we 
obtain for t = t : 



8nGp ^ r 2 2r 2H Q 

Hq ± R 2 Hq ± R'RHq ± h ± 



^0 = ^ + tttttt" + 7^772- - ^ = 1, (38) 



2 We do not provide any solution to the cosmological constant problem, namely, why A = 0. 
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and 

v 1 47rGp (r,t ) 1 r 1 H 0r {r,t ) 1 r 2 

glr ' r ° j 3 + 3H$ L (r, t ) + 3 RR>H* ± 3 # ±(r, t ) 3 R*H^ ' ^ 



If in (I39j) we have 



2 

flor^u. + ^ > 4vrGp + # 2 ± + (40) 



then the deceleration parameter g can be negative and cause t/ie universe to accel- 
erate without a cosmological constant or dark energy. If we solve for H 0r /H Q ^ from 
Eq.(j2HJ) and substitute into ()39|). we obtain 

8(r " tB) = 5-i(5i> (41) 

which gives a negative g for 

r 2 /i? 2 ff 2 ± > 3/5, (42) 



and we can satisfy the condition (|58|). We see that for r = corresponding to 
f 2 {r) = 1 and a spatially flat universe q(r,to) = 1/2. For p = A = 0, Eqs.(J7|) and 
(|40p lead to R > as the condition for an accelerating expansion of the universe. 

The conditions (f3*Kj) and (pfUJ) require that f 2 (r) > 1 for the inhomogeneous en- 
hancement, corresponding to an under-dense void. If we choose instead to describe 
the inhomogeneity by the spatially flat solution f 2 (r) = 1, then we cannot simulta- 
neously satisfy = 1 and q < 0. However, we will investigate in the next Section, 
whether the condition (J4~U|) can indeed be satisfied for the particular LTB model 
that we have adopted. We shall find that a more general inhomogeneous solution is 
needed to establish whether such a solution can explain the accelerating expansion 
of the universe without a negative pressure fluid or a cosmological constant. 

Several authors have investigated the robustness of the distance of the FLRW 
model determination of the magnitude red shift relation and the determination of the 
luminosity distance D L as a function of the parameters z, H , fl M and |5U|I5T| I2T]. 
In particular, an analysis of the magnitude red shift relation by Celerier [21] reveals 
that there exists a significant degeneracy and uncertainty in the determination of 
these parameters when analyzed using a LTB model to describe a local inhomogene- 
ity. 



5 The Raychaudhuri Equation 

We must now investigate whether the conditions Eqs. (j40|) and f!42|) can be satisfied 
by our LTB model. We will find that this is not the case and that we have to 
consider more general inhomogeneous solutions. To this end, we shall investigate 
the consequences of the Raychaudhuri equation PBICS], which holds in general for 
any cosmological solution based on Einstein's gravitational field equations. 
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We shall consider a congruence of curves with a time-like unit vector with 
g^VyVy = 1 and dV^/ds = V u V tl V 1 ' is the acceleration of the flow lines. The metric 
tensor W is given by 

h\ = 5\ + V»V v (43) 

and describes the metric that projects a vector into its components in the subspace 
of the vector tangent space that is orthogonal to V. 
We define the vorticity tensor 

fV = V^j, (44) 

and the shear tensor 

<V = V _ g^/UO ( 45 ) 

where 

V = V(^). (46) 
The covariant derivative of V can be expressed as 

= ^ + ^ + ^0 - ^7^- ( 47 ) 
The volume expansion is given by 

6 = hTQ^ = V^. (48) 
The Raychaudhuri equation is 

^ = -R V»V V + 2uj 2 - 2g 2 - -6 2 , (49) 
as 3 

where we have adopted geodesic world lines and 

uj 2 = uj^uj^ > 0, a 2 = g^g^ > 0. (50) 
From Einstein's field equations, we have for p = A = 0: 

f- = 2u 2 - 2a 2 - -9 2 - AnGp. (51) 
as 3 

The deceleration parameter q is defined by 

(Ue/ds + 6 2 ) _ 6(g 2 -u 2 )+4nGp 
q = W 2 ~~ 02 ' ^ 52 J 

If the vorticity uj = 0, then the deceleration parameter q has to be positive for 
positive p. 

Let us now consider the inhomogeneous metric (J2J) and the generalized Friedmann 
equations © and (J7J). The metric is chosen to satisfy the synchronous gauge and 
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comoving conditions goo = 1 and g 0i = (i = 1,2,3) corresponding to = 
(1,0,0,0). We have 

V M V" = d^V + T\ a V a . (53) 

We obtain for the metric (j2J: 

= rv = rv (54) 

This in turn becomes (see Appendix A in |T|): 

6 = ^ + 2 f = Hr + 2H± ' (55) 



For our synchronous comoving coordinates we have 

VVpO = V^d^O = 6. (56) 



de 

ds 



From (|53j) we get 



(36 + 9 2 ) \ 3 7v+ 6 r + 9 - 3H r ~ qh a ; 
9 - = p ( 57 ) 

We may now conclude that if 

m 2 r + m\ > 3^ + 6^ + e\ (58) 

then q could be negative leading to a local acceleration of the universe. However, 
for our synchronous comoving frame, we have 

= V W V A = d v V» - dpVp = 0, (59) 

so that the vorticity for this chosen gauge is zero. Therefore, we must deduce 
that the condition (|5Sj) cannot be satisfied for the positive energy condition p > 0. 
Moreover, the two conditions ()40j) and (|42j) cannot be satisfied for the synchronous 
and comoving metric (J2J). However, in the following section, we will show that 
we must carry out a spatial volume averaging of the expansion parameter 9 and its 
time evolution to arrive at a physically viable description of the local inhomogeneous 
late-time universe and an accelerating expansion of the universe. 



6 Spatially Averaged Cosmological Domains 

For our inhomogeneous model, it is necessary to perform a spatial average (|37|) of 
physical quantities, due to their observer, location dependence. If we postulate that 
the universe satisfies the cosmological principle and spacetime is described by the 
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homogeneous and isotropic FLRW metric, then we can make global predictions about 
the universe within our Hubble horizon. However, in general we cannot assume 
that the universe is homogeneous and isotropic, so we determine the structure of 
spacetime from the cosmological data available in our past light cone. By adopting 
a "smoothing" or averaging of the granular structure within a domain D in the past 
light cone, we can infer the properties of spacetime using this averaging process 
outside our past light cone (321 ESI El| • 

Although the gravitational potential $ is weak, $ < 1, we have to treat any 
back-reaction effect due to inhomogeneities in a relativistic way accounting for the 
non-linearity of Einstein's field equations. Sizable post-Newtonian effects due to an 
acceleration of the universe will appear in the Friedmann equations, describing an 
inhomogeneous universe after the smoothing or averaging process has been carried 
out. After a choice of time slicing is made to carry out the volume spatial averaging 
over a domain T>, the question whether an observer sees the universe accelerate or 
not depends on the data set the observer uses and the size of the scale over which 
the averaging process is invoked. Our model for the accelerating universe is based 
on a local void which plays the role of a cosmological constant. An observer living 
near or inside the void sees the local region of spacetime expand faster than the 
outer region of the void, despite the fact that both regions are decelerating. Thus, 
even though individual local fluid elements experience a deceleration, the averaged 
domain including the void and its outer regions will appear to the observer as an 
accelerating universe. 

An open question regarding the averaging procedure is whether it is covariant. 
The questions of covariance, gauge invariance and the uniqueness of the averaging 
procedure remain to be solved. The Einstein field equations are non-linear, so that 
the average of the tensor density of matter does not commute with the average 
of the time derivatives of the metric tensor occurring in the Einstein tensor. The 
average of the expansion of the universe derived from the luminosity distance must be 
consistently matched with the results obtained from the volume expansion parameter 
9 and the average of the density of matter. 

Let us define a more specific domain averaged expression for a scalar quantity [321 



is the volume of the simply-connected domain, D, in a hypersurface. We can define 
an effective scale-factor for our spatially averaged spherically symmetric inhomoge- 
neous enhancement (void): 




(60) 



where 




(61) 




(62) 
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where Vw is the initial volume. 

The volume averaging of the scalar \1/ does not commute with its time evolu- 
tion [321 1331 EH: 

(*(r,t)) D - dt(*(r,t)) D = (*(r,t)) D (9(r,t)) D - (*(r,t)9(r,t)) D . (63) 
We can derive for the averaged (9) d the equation 

(0)d = 3^ = S(H ± ) D . (64) 

Setting ^ = 9 and substituting for the Raychaudhuri equation (JoTj) for an irrota- 
tional matter dominated late-time model with uj^ u = 0, d9/ds = d9/dt and A = 0, 
we obtain the expression 

{R)d 

where (P)d is a function of (9)d, the average Hubble expansion parameters (H±) D 
and {H r ) D and the average shear (er) D . 

For inhomogeneous cosmology, the smoothing due to averaging of the Einstein 
field equations does not commute with the time evolution of the non-linear field 
equations. This leads to extra contributions in the effective, averaged Einstein field 
equations, which do not satisfy the usual energy conditions even though they are 
satisfied by the original energy-momentum tensor. It is the lack of commutativity 
of the time evolution of the expansion of the universe in a local patch inside our 
Hubble horizon that circumvents the no-go theorem based on the local Raychaudhuri 
equation. 

If we can satisfy the condition 

(P) D > AnG(p) D , (66) 

then the averaged deceleration parameter (q) d can be negative corresponding to 
the acceleration of a local patch of the universe. If we substitute the value of (R) d 
from the volume averaged equation (JJJ) for p = A = 0, then the exact LTB solution 
for an irrotational inhomogeneous, under-dense void solution can describe the local 
acceleration of the late-time universe and fulfill the constraint equation ()66j) . The 
spatial domain average of q in Eq.(|4ip gives 

(<,)o=\-U(^A) . (67) 



/ J2. 



2 

This yields the averaged condition for an accelerating universe: 

r 



\Rmz 



> 3/5. (68) 



'0_L / D 

This spatial domain averaged condition for our LTB void solution can be satisfied 
in our synchronous comoving gauge with zero vorticity for (p) D > 0. 
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7 Conclusions 



We have modelled the late-time inhomogeneous and non-linear regime for z < 10 by 
an exact spherically symmetric, inhomogeneous solution of Einstein's field equations. 
The solution approaches an Einstein-de Sitter matter dominated solution for z > 10, 
which describes a homogeneous and isotropic FLRW background spacetime with 
small inhomogeneities 5p/p ~ 10~ 5 near the surface of last scattering in agreement 
with the WMAP3 data jlj. An investigation of the Raychaudhuri equation for the 
expansion rate of the universe, shows that for the synchronous and comoving gauge 
condition satisfied by the inhomogeneous spherically symmetric metric, the vorticity 
uj is zero and the deceleration parameter q is positive or zero when the positive energy 
condition is satisfied. However, a volume averaging of physical quantities is required 
for our inhomogeneous universe, and due to the lack of commutativity of the time 
evolution of the average of the scalar expansion 6, we find that it is possible to have 
a negative averaged deceleration parameter (q) for a positive density of matter p and 
a zero vorticity uj = 0. Observational bounds on the magnitude of the quantity (Pd) 
in the condition (jfifij) must be obtained to establish whether a sufficiently negative 
q can be produced by the expanding inhomogeneous enhancement or void to allow 
an explanation of the accelerating universe without a negative pressure fluid or a 
cosmo logical constant. 

The model inhomogeneous solution is exact, so that we do not have to be con- 
cerned with the failure of perturbative backreaction calculations in the non-linear 
late-time universe regime. However, we have assumed a high degree of symmetry 
for our inhomogeneous solution, since it only allows for a one- dimensional inhomo- 
geneity in the radial direction, although the angular azimuthal expansion described 
by H± does play a significant role in determining the cosmological solution and the 
magnitude of the acceleration. Due to the need for approximate isotropy, we re- 
quire that the observer be situated not too far from the "center" of the spherically 
symmetric, inhomogeneous distribution of matter, so this could be interpreted as a 
"coincidence" problem. But this could be an artifact of the degree of symmetry as- 
sumed for the model solution. It would be interesting to investigate inhomogeneous 
cosmological solutions that possess less symmetry than the LTB models j29| . 
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